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Abstract 

The paper is concerned with thermostatistics of both D-dimensional Bose and 
Fermi ideal gases in a confining potential of type Ar n + Br~ n . The investigation 

•i-H . 

is performed in the framework of the semiclassical approximation. Some physical 
quantities for such systems are derived, like density of states, density profiles and 
number of particles. Bose-Einstein condensation (BEC) is discussed in the high and 
low temperature regimes. 
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1 Introduction 



In the recent past there has been increasing emphasis in quantum statistics of the 
ideal quantum gases in an arbitrary quantum potential since the recent experiment obser- 
vations. The latter concerns the Bose-Einstein condensation |lf and the Fermi quantum 
degeneracy H for dilute alkali-metal atoms in magnetic or magneto-optical traps at very 
low temperature. These results have received renewed attention and much of this atten- 
tion is paid to the study of ideal quantum gases confining in different external potentials. 
In particular, the Bose-Einstein condensation. This was examined in many occasions in 
the presence of typical external potentials, for example; harmonic potential [3—8], toroidal 
potential 0, double- well potential |10| and in the presence of an impurity [11]. More re- 



cently, an interesting results have been obtained by Salasnich. Ideal quantum gases in 
.D-dimensional space and confining in power-law potential fl2| have been studied. Among 



the results derived in this paper, we note the condition of the Bose-Einstein condensation 
(BEC). Indeed, it is shown that BEC can set up if and only if y + — > 1. 

In this paper we wish to extend the Salasnich analysis fl2[] to cover more general 



results. For this task, we investigate the ideal quantum gases trapped in the potential 

U(r) = Ar n + Br~ n (1) 

with A, B strictly positives and n is a non vanishing integer. For this type of poten- 
tial, we give the thermostatistics properties which generalize the Salasnich ones. In the 
limit B — ► 0, we recover the Salasnich results. We prove that BEC occurs in the high 
temperature regime. 

This paper is organised as follows. In section 2, we give the quantum distribution func- 
tions in the framework of the semiclassical approximation for both bosons and fermions 
confining in a generic potential. Section 3 is devoted to the generalization of the Salas- 
nich analysis when the particles, bosons and fermions, are trapped in the potential given 
by equation (1). We will see that in the limit B — ► 0, one can recover the results ob- 
tained by Salasnich in . Bosons and fermions are studied separately. For the fermionic 
gas evolving in the U(r) potential, we give the corresponding Fermi functions, the finite 



and zero temperature momentum distribution, and the number of particles. In a similar 
way, we define the Bose function corresponding to particles obeying the Bose-Einstein 
statistics, we compute the finite temperature non condensed momentum distribution and 
finally we give the relation satisfied by the Bose transition temperature below which the 
BEC takes place. Since, this relation can't be solved in the general case, we restrict the 
BEC analysis in the high and low temperature regimes. The last section is devoted to 
the conclusions and perspectives of the present work. 



2 Prelimenaries 

We shall briefly review the basic tools which will be useful to investigate the thermal 
properties of both Bose and Fermi ideal gases in a generic confining external potential. 



2.1 Quantum distribution functions 

We start with a confined quantum gases (identical bosons or fermions). The average 
number of occupation for the k-mode is [[DJ 

N * = e/»(«*-M) + l' (2) 



for the fermions and 



for bosons. In equations (2) and (3), fi is the chemical potential, are the energy of 
particles, and (3 = -j^p the reciprocal temperature. Then, the total average occupation 
number iV of the particles is given by 

# = (4) 

k 

which leads to fix the parameter \i. For more details about the physical meaning of the 



chemical potential \i see ref.[|12 . 
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In conclusion, equations (2) and (3) will be written in the semiclassical approximation 
in the next subsection and will be constituted a starting point of our study for the bosons 
and fermions, respectively. 



2.2 Semiclassical approximation 

In the semiclassical approximation (the thermodynamical limit), the energy spectrum 
may be considered as a continuum. This situation occurs when the number of particles 
is large and the energy level spacing is small. Thus, the quantum distribution functions 
can be replaced by the so-called phase-space distribution JT^ 

^P)= e ^ (r ,p)-,) ±r (5) 

where r = (r\,r2, , tq) and p = (pi,P2, ,Pd) are the position and momentum of D- 

dimensional system under consideration. In equation (5), plus and minus refer to fermions 
and bosons respectively. To obtain the total distribution, we compute the integral over 
all 2D-dimensional phase-space. Then, we get 

r d D rd D p , . 

iV = /w" (r ' p) ' (6) 

which can be written as follows 

N = J d D r n(r) = / d D p n(p), (7) 
where the spacial and momentum distributions are given, respectively, by 

n(r) = /(4^" (r ' p) ' (8) 



and 



n(p) = /(2^" (r ' P) - (9) 



Introducing the density of states 



p(e)= /f^F^-^'P^' (10) 
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with 5(x) is the usual Dirac function. The average total occupation number of particles 
can be written as 

The phase-distribution, in the case of fermions at zero temperature (i.e where the chemical 
potential /i coincides with Fermi energy Ep), takes the simplest form 

n(r,p)=e(E F -e(r,p)), (12) 

where Q(x) is the well-known Heaviside function. The case of bosons was discussed in 
great detail in |l2j (see also the references therein). When, the external potential is taken 
in to account, the classical single-particle energy is defined as 

^,p) = f^(r), (13) 
from which one can prove that the semi-classical density function is given by 

^-^'fk)I dDr{e - u{r)) ^- (14) 

where T(n) is the Gamma function. It is interesting to note that equation (14) is valid for 
all external potentials. We now close this section related to the main tools which will be 
used in the following sections investigating the thermostatistics of bosons and fermions 
trapped in the potential given by eq.(l). 

3 Ideal gases in the Ar n + Br~ n potential 

In many experiments with alkali-metals atoms, the external trap can be modelled by a 
harmonic potential |TJ]]. The effects of adiabatic change in the trap can be represented 



by power-law potentials |L2[]. In general, to find the momentum distribution, the Fermi 



temperature and the Bose temperature, it is necessary to specify the external potential. 
3.1 Fermionic gas 

To discuss the statistical properties of a fermionic gas embedded in the potential U (r) , let 
us start by defining the Fermi function. As we will see, this definition generalizes the one 



given in [1^] and allows the computation of physical quantities, of interest in the study 
of the system under consideration, as, for instance, the density of states and number of 
particles. 

We define the Fermi function as follows 

f™W = rM / dyy*- 1 (15) 
1 [n) Jo 1 + ze ay oy 

where a and b are two arbitrary positive constants. In the particular case b = 0, we get 
after rescaling the variable 

I r°° 7p~ y 
Uz) = — / dyy^ 1 -^-— , (16) 
1 in) Jo 1 + ze v 



where f n (z) is the Fermi function definition [T!| for a fermionic system in the U(r) = Ar n 



potential. We note that, the Fermi function can be expanded, for \z\ < 1, as 

fi B \z) = YTV £ * K n (2yfib), (17) 

in terms of the modified Bessel function K n (2y/ab) [l5j 

1 f, POO 

K n (2Vd>) = -(-)? / dyy n -\- a ^-\ 



o 



The density of states of a fermionic quantum gas in the potential eq.(l) can be calcu- 
lated from equation (14). So, we have to compute the integral of type 



P B ( e ) = (t^) T T^TJT / d ° r ( e - Ar " - Br- n )^. (19) 



m d 1 

2 F( 2 , .. 

In this order, we put y = r n and solve the following equation 

e - Ay - By' 1 = 0. (20) 
The latter equation has two positive solutions 



e / . € . „ B e / , 6 s „ B . . 

"=2l + V ( 2l )2 -^ tt =2l-V ( 53»-^ (21) 
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when e > 2 V ' AB is satisfied since the radius of the sphere r should be a real positive. 

D 

Using the relation d D r = giving the D-dimensional unit sphere and the equations 

J- ( 2 +1) 

(21), the integral (19) can be written as 

D . D-2 

R . . . Tfl . D DTI 2 A 2 fU D D , . D-2 , , D-2 . . 

P e) = ( — t)~ ,m 7n r / dy y-~-(a + y) — (u-y) — , (22) 

F w V 27r^ ; n r(f)r(f + i) y« v ; v ; K J 

which can be solved using the integral representations of the hypergeometric functions. 
Finally, we obtain 

r> . . . m .d D d-2 . n-2 r(£-g + l) . D D D D u s . , 

' (£) = < va» T n W T [r ( f + V " 2 ^ " 2 + 1; 7T ; (23) 



F 1 (a,ft 7 ,^)= „, ' „ r^/^d-y) 7 """^!-^)"", (24) 



We recall that the hypergeometric function 2 -Fi(o!,/3; 7; z) is defined by JT5 

! / 

B(/3,7-/3) Jo 
where _Re7 > 0, Re/3 > and 

^ a '^ = l{T+p) =ma) ' (25) 

is the so-called beta function or Euler's integral of the first kind, which admits the following 
integral representation 

B{a,/3)= [ dtt^il-tf- 1 . (26) 
J 

It is interesting to remark that the following relation is satisfied 



, , xA , D D D D u, ,e,D , T(%)T(£ + l)r(^) 



(27) 

Using the latter relation, the density of states, in the limit B — > 0, is given by 

, s , Tfl . D , \ , D r(— + 1) D ,D 11 

»< £ > = <5? )T <3>* r ( f + 1 V(f + f) ^~ (28) 



which coincides with the result obtained in reference 

Using the above tools, other interesting properties of the present system (fermions) 
may be determined in a simple way: 
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i- The finite temperature momentum distribution 



ii- The zero temperature momentum distribution 

where we request that p < [2m(e — 2\/AB)]5 for the right term makes sense. 

iii- Number of particles 

The number of particles is given by the following integral 

N = * / jgP. q {Ef - ^) f F ~ & + [( ^ " - -]¥ 
T(f + 1) 7 (2tt^) d 1 F 2m J [ 2A ^ [[ 2A } A 1 1 ' 

which can be reorganized as 



(31) 



~2 _ „2 



1 DtT~ r[2m(e-2VAB)]? E F - w E F -2-^ B,i,d 

N -(wm^wl d " p -'i-2?* + i<--2?*y-AW- 

(32) 

A calculation more or less complicated, gives the following expression for the number of 
particles 

N = M /( 7 ) B(%, [A ( 7 ) 2 F 1 (+1) - /f 1 ^) 2 F 1 (-1)], (33) 



where 



2 ' 2 

/( 7 ) = ( 7 2 -l) — [ 7 -(7 2 -l) 5 ]-"-, 
A(7) = 7-(7 2 -l)i 

/r 1 (7)= ! 



E F - " 

7 = 



,2 



2A ' 

, . . D D D 2( 7 2 -l)5 

2 F 1 (±1) = 2 FA + — =f 1; — ; D, '—), 

2 u ; 2 u n 2 ^ ' 2 ' ' 7 _ ( 7 2_ i)^' 



and £(§, f ) is given by eqs.(25,26). 



The number of particles is expressed in terms of the hypergeometric and beta functions. 
It is clear that this relation is complicated and we believe that a more detailed study of 
it can generate more interesting results. We will return to this matter in a subsequent 
paper. 



To close this subsection, we note that the results derived in ref.p2| related to the 
fermionic g be reproduced just by taking the limit B — > 0. 

3.2 Bosonic gas 

We now examine the condensation of a system of bosons confining in the potential U (r) 
(eq.(l)). We start by defining the Bose function as follows 

1 /-oo 7p -ay-by~ 1 
g (B)^ = " / dyy n-l_^ ( 34) 

yn v ; r(n) Jo y y i - ze-w-ty- 1 v ; 

For \z\ < 1, we have 

9n B \z) = E K n (2jV^b), (35) 



where K n (2y ab) is given by eq.(18). For B — > 0, we show that 

9n{z) = —- / dyy n - 1 -^— , (36) 
i[n) Jo 1 + ze y 

which is nothing but the definition of the Bose function given in corresponding to the 
bosons in the potential U(r) = Ar n . 

Using the definition given by eq.(34), one obtain the following results: 
i- The finite temperature non condensed momentum distribution 

2 



n(p) " {2jRK)° r(f + 1) 9r {e ' 1 

ii- The Bose transition temperature 



To obtain the Bose-temperature, we have to calculate the number of particles 

N = J d D pn(p). (38) 
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As in the case of ordinary bosons (without interaction), we define the Bose-temperature 
by taking the chemical potential \x = 0. By introducing \i = in eq.(37), we get 

g*j r(g + i) 
(2^h) D [r(§ + i)p 

With the help of eq.(35), it is easy to show that 



n = r ;;j / ^-^(^H ( 3 9) 



N = 2( — 135 



R D ,2m N D D7r^ rf^ + l) ~ K R (2jPVAB) 

\ 7^ — / \ TT \ 71 / \ 11 



77 E £ • (40) 

1 _•_ 1 7 2 



A 7 v /3 ; {2^h)° r(f + i)r(f + . =1 j 2 

The latter equation gives a relation between the number of particles (or density) and 
Bose-temperature below which we obtain a Bose-Einstein condensation. Due to the com- 
plicated sum appearing in (40), the Bose temperature can't be calculated, in general, in a 
simple way. However, in the some particular situations, the solutions of this equation can 
be found. Indeed, in the high and low temperature regimes (which are two interesting sit- 
uations from an experimental point of view) some information about BEC can be derived: 



High temperature: This case corresponds to (3 — ► 0. In this limit the functions Kd {2j[3\J AB) 



take the form [ |i~6|j 

,D 



Kn(2j(3VAB) « - r(-) (2j (3VAB)--. (41) 
2 n 



N - [ 2tf ] r(f) ^#+f ' (42) 



Therefore, eq.(40) implies 



This result is similar to the one found in ref. ||12|| . From this relation we conclude that just 
restricting ourselves to high temperature limit, we find the Salasnich analysis related to 
the BEC. Indeed, it is shown that the latter takes place when the condition + — > 1 is 
satisfied. Strictly speaking, we obtain the condensed fraction as follows 

Nn , T . D ,D , 

^ = l-(^, (43) 

where Nq is the number of particles occupying the single-particle ground-state of the 
system when the temperature is below Tg. 



Low temperature: Equivalently to f3 — > oo, then we have |15 



Ko(2j(3VAB) « 0. (44) 

n 
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In that case, we get a number of particles N — > corresponding to a temperature 
T — ► 0. This result is compatible with the literature |jl3[ . 



4 Conclusion 

In this work, we concentrated on the influence of confining potential U (r) = Ar n + Br~ n 
on the main statistical properties of the ideal quantum gases (bosons or fermions) in D- 
dimensions. We discussed the derivation of the density of states, spacial and momentum 
distributions in the thermodynamical limit. For fermions, we have calculated the Fermi 
energy and for bosons, the phenomenon of Bose-Einstein condensation is discussed in 
terms of the reciprocal temperature ((3). Two situations were considerd ((3 — > 0) and 
{(3 — > oo) corresponding to high and low temperature, respectively. In the limit B — ► 



0, our results reproduce the Salasnich ones ref.[12| concerning the confining power-law 
potential. 

An important result of the present work concerns the Bose-Einstein condensation, in 
the high temperature domain, which occurs when ^ + — > 1, where D is the space di- 
mension and n the exponent appearing in the expression of U(r). Another result concerns 
to the low temperature case where the number of particles vanishes. 

To finalize this paper, it should be noted that the present investigation of the ideal 
quantum gases embedded in the potential U(r) (eq.(l)) in D-dimensions, constitutes now 
a very interesting topic to learn more about Bose-Einstein condensation. In fact, as we 
have indicated above, the BEC should be investigated in more detail in general case 
solving equation (40). We believe that more important informations can be obtained in 
this case. We will return on this subject in a forthcoming work |jT^ . 
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